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The key characteristic of an optical mode in a microcavity is its quality factor describing the 
optical losses. The numerical computation of this quantity can be very demanding for present-day 
devices. Here we show for a certain class of whispering-gallery cavities that the quality factor is 
related to dynamical tunneling, a phenomenon studied in the field of quantum chaos. We extend 
a recently developed approach for determining dynamical tunneling rates to open cavities. This 
allows us to derive an analytical formula for the quality factor which is in very good agreement with 
full solutions of Maxwell's equations. 

PACS numbers: 42.55.Sa, 42.60.Da, 05.45.Mt 



I. INTRODUCTION 

Optical microcavities in which photons can be con- 
fined in three spatial dimensions are a subject of inten- 
sive research as they are relevant for applications, such 
as ultralow-threshold lasers [2j , single-photon emitters 
[H, El or correlated photon-pair emitters [|| . Especially 
whispering-galler y c avities such as microdisks p, 0i @|j 
microspheres [9|, [lOj, and microtoroids [ll| have been 
investigated as they can trap photons for a long time 
near the boundary by total internal reflection. The cor- 
responding whispering-gallery modes have a very high 
quality factor Q, which makes these cavities a candi- 
date for the above-mentioned devices. While the mi- 
crodisk emits the photons isotropically, cavities with de- 
formed surfaces may additionally lead to directed emis- 
sion [n m q m m, m m m, m m m m. 

A particularly interesting geometry is the annular cav- 
ity [li], [25|, H|| - a microdisk with a circular-shaped in- 
clusion. A non-concentric (air) hole as inclusion allows 
for unidirectional emission and high quality factors simul- 
taneously [13] which for most applications are desirable. 

In this paper we connect the quality factors Q of opti- 
cal microcavities to the concept of dynamical tunneling 
f27l ]. We provide an explicit prediction for the quality 
factors of whispering-gallery modes in the annular micro- 
cavity; see Fig. [T] Microcavities typically have a mixed 
phase space where regions of regular and chaotic mo- 
tion coexist. Dynamical tunneling occurs between these 
dynamically separated phase-space regions. While for 
one-dimensional systems the tunneling process through 
an energetic barrier is well understood, e.g., by means 
of WKB theory [28|, dynamical tunneling is a subject 
of intensive research experimentall y I29L 130 . 13111 as well 
as theoreticallyjp [H, M, M,M,M,M, M, IM EE 
EH, IHj EE El; l45j. Using a fictitious integrable system 



[4l|, recently an approach has been developed [44|, |45( 
which successfully predicts dynamical tunneling rates for 
quantum maps and billiards. In this paper we extend 
this approach to open optical microcavities, in particu- 
lar to the annular cavity, to predict quality factors Q of 
whispering-gallery modes. 

This paper is organized as follows. In Sec.|TT]we intro- 
duce the annular cavity. In Sec. IIIII we define the quality 
factors and derive the dynamical-tunneling contribution, 
which is then compared with numerical results. A sum- 
mary is given in Sec. IIVI 
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FIG. 1: (Color online) Inverse quality factors 1/Q for the 
annular microcavity with refractive index n e ft = 2.0. Shown is 
the theoretical prediction (solid line) , which is the sum of the 
direct tunneling contribution (dotted line) and the dynamical 
tunneling contribution (dashed line, Eq. (|15[l l. and numerical 
data (dots) for angular quantum number m = 19 and radial 
quantum number n = 1 vs the hole position d. The insets 
show the resonant state at d = 0.2475 and d = 0.33. 
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II. ANNULAR MICROCAVITIES 

Optical cavities are described by Maxwell's equations 
which in the case of quasi-two-dimensional microdisks 
reduce to a two-dimensional scalar mode equation [46| 



V 2 ^ = n 2 (x,y)k 2 if} 



(1) 



with (effective) index of refraction n(x,y), wave number 
k = ijj/c, frequency u>, and the speed of light in vac- 
uum c. The mode equation (fTJ) is valid for both trans- 
verse magnetic (TM) and transverse electric (TE) polar- 
ization. We focus on TM polarization with the electric 
field E(x, y, t) = (0, 0, ip(x, y)e~ luJt ) perpendicular to the 
cavity plane. The wave function ip and its normal deriva- 
tive are continuous across the boundary of the cavity. At 
infinity, only outgoing-wave components are allowed. 

The mode equation {!]) with the above-mentioned 
boundary conditions has analytical solutions only for spe- 
cial geometries such as the circle (see the Appendix) 
or several concentric circles. General geometries re- 
quire numerical schemes such as the boundary element 
method [47| which we use in this paper. In this approach 
the two-dimensional partial differential equation ([T]) is 
rewritten as a one-dimensional integral equation involv- 
ing Green's functions. Outgoing-wave conditions can 
be easily fulfilled by using the outgoing solution for the 
Green's functions. The boundary element method turns 
out to be very efficient even for computing highly excited 
modes and their quality factors. 

A particularly suited example to study the influence 
of a mixed phase space onto the quality factors Q is the 
annular cavity. Its geometry is given by the radius R of 
the large disk, the radius R2 of the small disk and the 
minimal distance between the two disks d, see Fig. |2£a). 
Without loss of generality we choose R = 1. Under 
the variation in the two parameters d and R2 the dy- 
namics inside the cavity changes drastically from com- 
pletely regular behavior when the two disks are concen- 
tric (d = R — R2) to mixed regular-chaotic behavior for 
the general eccentric case. This is clearly visible for dif- 
ferent trajectories in the Poincare section, see Fig. [2Jb). 
The Poincare section is a two-dimensional phase-space 
representation. Whenever the trajectory hits the cav- 
ity's boundary, its position s (arclength coordinate along 
the circumference) and tangential momentum p = sin% 
(the angle of reflection x is measured from the surface 
normal) is recorded. The large disk has an effective re- 
fractive index n e g while inside the small disk and outside 
of the cavity the refractive index is unity. For the visu- 
alization of the ray dynamics we used an annular cavity 
with hard wall boundary conditions at the outer disk, 
neglecting ray-splitting effects. 

The annular cavity has been studied extensivel y in 



26], 
24 j. 



the context of quantum chaos [25l ] . optomcchanics 
avoided resonance crossings and resonant tunneling 
For applications it is of high interest as it allows for uni- 
directional light emission from high-Q modes which has 
been predicted by the authors in Ref. [171 ] and has been 



confirmed in recent experiments [2l|. The closed system 
with perfectly reflecting walls, i.e., the annular billiard, 
is a paradigm for dynamical tunneling (3f| HH . 



III. QUALITY FACTORS 

The quality factor Q of a mode in an open cavity is re- 
lated to the corresponding resonance with complex wave 
number k = Rc(fc) + ilm(fc) via 



Q = - 



Re(fc) 
2Im(fc) ' 



(2) 



In the annular cavity the quality factor Q of a regular 
mode has two contributions which are assumed to be ad- 
ditive 
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/dir 



Qdyn 



(3) 



Here, Qdir accounts for the direct coupling of the regular 
mode to the continuum, as in the case of the circular 
cavity, see the Appendix. Note that for this contribution 
the mixed phase-space structure induced by the small 
disk is irrelevant. The second contribution, Qdyn, is given 
by dynamical tunneling from the regular mode to the 
chaotic sea, which for \p\ < p c is strongly coupled to the 
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FIG. 2: (Color online) (a) Annular cavity with R — 1, 
R2 — 0.22, and d — 0.33. (b) Corresponding Poincare sec- 
tion (s,p = sinx) for n e ff = 2.0, where s is the arclength 
along the boundary and x ls the angle of reflection. Between 
the critical lines (dashed) with |p c | = l/n e fl is the leaky re- 
gion, where the condition for total internal reflection is not 
fulfilled, (c) Circular cavity and (d) its completely regular 
phase space. 
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continuum, see Fig. [2f b) . Here we assume that there are 
no further phase-space structures within the chaotic sea 
that affect the quality factor. A priori it is not obvious, 
which of these contributions will dominate. 



A. Dynamical tunneling contribution 

We now want to derive a prediction for the dynamical 
tunneling contribution Qdyn of regular modes in optical 
microcavities. After presenting the general approach we 
will apply it to the annular microcavity. 

1. General approach 

We first review a quantum mechanical approach for 
determining dynamical tunneling rates using a fictitious 
integrable system [44|, |45| ; the relation to the quality fac- 
tor Qdyn will be given below. The tunneling rate 7 of a 
regular state to the chaotic sea is described by Fermi's 
golden rule (using units h = 2M = 1) 

7 = 27r(M) 2 /9 ch (4) 

where p c h denotes the density of chaotic states and (\v\) 2 
is the averaged squared matrix element between the con- 
sidered regular state and the chaotic states of similar 
energy. According to the Weyl formula for closed two- 
dimensional billiards the density of chaotic states is given 
by Pch ~ Adi/(47r), where A c ^ is the area of the billiard 
times the fraction of the chaotic phase-space volume. The 
eigenmodes of a system with a mixed phase space are 
mainly regular or chaotic, i.e., concentrated on a torus 
inside the regular region or spread out over the chaotic 
component. To calculate the coupling matrix elements v 
these so-called regular and chaotic eigenmodes cannot be 
used as they have small, but still too large, admixtures of 
the other type of modes compared to the tunneling rate. 
Instead, we determine v by introducing a fictitious inte- 
grable system H reg as it was first suggested for dynamical 
tunneling in Ref . [4l| . H Teg has to be chosen such that 
its classical dynamics resembles the regular dynamics of 
the mixed system as closely as possible and extends it 
to phase space regions where H has a chaotic sea. The 
eigenstates V>reg of H leg are localized in the regular region 
of H and decay into the chaotic sea of H . With chaotic 
states V'ch, which live in the chaotic region of phase space, 
the coupling matrix element is given as [3 EH, |48| 

v = i>*h( x i V){H - H rcg )ip rcg (x, y) dxdy. (5) 
Jr 2 

Note that this equation is applicable for general systems 
but the determination of a sufficiently accurate H Teg is a 
difficult task. As previously mentioned we assume that 
there are no further phase-space structures within the 
chaotic sea that affect the tunneling rates 7. This ap- 
proach was previously used to predict tunneling rates for 
quantum maps [H| and closed billiard systems [45| . 



The described approach can be extended to open cav- 
ities in the following way: (i) as a fictitious integrable 
system H rcg we choose a cavity such that it resembles 
the regular dynamics of H. The quantum system has 
resonance states ipreg- (h) As a model for the chaotic res- 
onances ipeh a random wave model will be used, which in 
addition fulfills the relevant cavity boundary conditions, 
(iii) The tunneling rate 7 determines the quality factor 
Qdyn by 

_ 21^^ 

7 

where we used the quantum mechanical relation between 
energy and momentum Eq — i-^ — p 2 , with p = n c ffk 
in a refractive medium such that 7 = — 2n 2 ff Im(fc 2 ) = 
-4n 2 ff Re(A;)Im(fc). 

2. Application to the annular cavity 

To evaluate Eq. ([5]) we have to find an appropriate 
regular system H rcg . For the annular cavity a natural 
choice is given by the circular cavity as it correctly repro- 
duces the regular whispering-gallery motion of the annu- 
lar cavity and extends it into the chaotic region of phase 
space. The whispering-gallery modes are labeled by the 
two quantum numbers m and n. Rewriting the mode 
equation ([T]) as an eigenvalue equation Hip — k 2 ip the 
Hamiltonian of the annular cavity can be introduced as 

H = -V 2 + [1 -n{x,y) 2 ]k 2 (7) 

where the refractive index n(x, y) is n e ff inside and 1 out- 
side the cavity and in the disk of radius i?2- As the reg- 
ular system H reg we choose the circular cavity 

ffreg = -V 2 + [1 - n rcg (x, y) 2 ]k 2 (8) 

where the refractive index n reg (x,y) is n e s inside and 
1 outside the circular cavity. Thus H and H rcg differ 
only inside the small disk of radius R2 with H — H Teg = 
(1 — n 2 s )k 2 and the integral in Eq. © reduces to an 
integral over the small disk. For the regular states inside 
the circular cavity we choose, as in the numerical studies, 
the even eigenmodes 

V>™g ( r : <P) = N mn J m (n cS k mn r) cos(rmp) (9) 

where k mn are the complex resonant wave numbers, ac- 
cording to the Appendix, and ?p™ g is normalized to one 
with the numerically determined normalization constant 

To model the chaotic modes ip c h within the small disk 
we employ a random wave description (49| . which has 
been extended to systems with a mixed phase space [5(| ■ 
While this model accurately describes the random behav- 
ior in a medium with constant refractive index, it cannot 
account for the change in refractive index at the border of 
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the small disk at p = Ri- We extend a boundary-adapted 
random wave model [5l| to account for this boundary 
condition. This is essential especially for low Re(£;) as 
then all chaotic modes decay inside the small disk, which 
cannot be reproduced by the usual random wave model. 
Therefore we construct the chaotic states ip^ as a random 
superposition of modes of a circular cavity of radius R2 
with refractive index 1 which is surrounded by a medium 
with refractive index n g. As follows from Eq. (fl"9|) these 
modes, at a fixed complex wave number k mn and within 
the small disk, are 



tpi(p, ■&) = AiJi(k mn p) cos(W), p< R 2 



(10) 



The chaotic states V'cii are then constructed by a random 
superposition of these modes 



V>ch(p, ■&) = —j=F= X! fli ^(P > 



(11) 



1 = 1 



Here the coefficients ai are Gaussian random variables 
with mean zero and (aiak) = Si^- The random waves, 
constructed in such a way, fulfill the normalization con- 
dition { IV'ch 1 2 ) = 1/A;h required for the annular cavity. 

Using the fictitious regular system and the random 
wave model for the chaotic states we obtain an integral 
over the small disk for the coupling matrix element 

R 2 i-T 



(12) 

For the tunneling rate this results in 
1 

Imn = ^N 2 mn {l - n 2 cS ) 2 \k 2 mn \ 2 J2 \!i\ 2 (13) 



1=1 



where 
h = 



p dpdtf tp* (p, •&) J m (n eS k mn r) cos{mtp) 



with r = r(p, ??) geometrically related to ip = <p(p, 1?), see 
Figs. [2^a) and (c). With Eq. ^ we finally obtain the 
dynamical tunneling contribution to the quality factor 

r>mn _ off /1 rs 

^ a Ni n (i-n 2 cS r\k 2 mn \j:Zi\ii\ 2 1 j 

of whispering-gallery modes in the annular microcavity 
for each quantum number m and n. 



B. Results 

Now we compare our theoretical prediction for the 
quality factor, Eq. ([3]), for the annular microcavity with 
numerical data, obtained using the boundary element 
method [13] • Figure [T] shows the inverse quality factors 



at fixed quantum numbers m — 19 and n — 1 under vari- 
ation in the distance d between the small and the large 
disk. The direct contribution l/Qdir, see the Appendix, 
is independent of the distance d. It is dominated by 
the dynamical tunneling contribution l/Qdym Eq. (|15p . 
which decreases exponentially with d, as expected from 
the increasing regular phase space region. We find excel- 
lent agreement of the prediction and the numerical data. 

As a further test we consider the quality factors for 
fixed radial quantum number n = 1 and increasing angu- 
lar quantum number m = 7, . . . , 21, comparing the theo- 
retical prediction with numerical results, see Fig. [3] We 
find that for small Re(fc) the direct tunneling from the 
whispering-gallery modes to the continuum is relevant 
while for large Re(fc) the dynamical-tunneling contribu- 
tion dominates. Here, our prediction again shows excel- 
lent agreement with the numerical data. Note that for 
other systems, such as the one considered in Ref. [42j], 
only Qdyn may be the relevant contribution. Also we 
point out, that our theory allows to determine quality fac- 
tors for large Re(fc), where numerical methods fail. The 
boundary element method cannot compute the quality 
factors of quantum numbers n — 1 and m > 21 reliably 
as the exponentially increasing quality factor requires an 
extremely fine spatial discretization. 



C. Additional phase space structures 

In the derivation of the dynamical tunneling contribu- 
tion 1/Qdyn to the quality factor we assumed that there 
are no further structures in the chaotic part of phase 
space, such as small regular islands and partial barriers. 
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FIG. 3: (Color online) Inverse quality factors 1/Q for the an- 
nular microcavity with n c ff = 2.0. Shown is the theoretical 
prediction (solid line) which is the sum of the direct tunneling 
contribution (dotted line) and the dynamical tunneling con- 
tribution (dashed line, Eq. (US}), and numerical data (dots) 
for m — 7, ... ,21 and n = 1 at d = 0.33. The insets exem- 
plarily show the resonant states of angular quantum number 
m = 7 (left) and m = 21 (right). 



If this assumption is not fulfilled, the tunneling rate 7 is 
modified and consequently the quality factor. 

This can be demonstrated when choosing d = 0.33, 
i?2 = 0.22, and increasing n g from 2.0 to 2.3. At 
n e ff = 2.0 no visible additional structures exist in the 
chaotic part of phase space above the critical line (see 
Fig. [2b)). At n c ff — 2.3 the critical line p c is shifted 
to smaller values and a period-three island chain is now 
above p c as can be seen in the inset in Fig. 0J These 
structures presumably cause the oscillations on top of the 
numerically determined quality factors that are visible in 
Fig. [U To support the conjecture, that the island chain 
is responsible for the oscillations, Figs.JSta) and[S]Jb) dis- 
play the incident Husimi functions [52| , representing the 
wave analog of the Poincare section (see also [HJ]), of 
the mode m = 14 (near the minimum of the oscillation) 
and m = 18 (near the maximum). In the former case 
the island chain is clearly a barrier for the mode. The 
mode cannot penetrate the leaky region so easily, which 
increases its quality factor. In the latter case the island 
chain seems not to have a strong influence on the mode. 
While the average behavior of the quality factors is still 
well predicted by Eq. (|15[) . these oscillations cannot be 
explained by our theory as it assumes a strong coupling 
of the chaotic modes to the continuum. Other situations 
where this coupling is weak are due to dynamical localiza- 
tion I54j| or to an additional tunnel barrier surrounding a 
cavity [2J| . We leave the interesting task of the prediction 
of tunneling rates through more complicated phase-space 
structures in the chaotic future challenge. 
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FIG. 5: (Color online) Husimi functions (gray scale and con- 
tour lines) of the modes with angular quantum number (a) 
m — 14 and (b) m = 18 superimposed onto part of the 
Poincare section of the classical phase space, where the criti- 
cal line p = p c is marked (dashed line). The radial quantum 
number is n = 1 and the refractive index is n c g — 2.3 as in 
Fig.H 



IV. SUMMARY 

We have presented a theory for the intrinsic optical 
losses of annular microcavities. It is assumed that the ray 
dynamical phase space is divided into regular regions and 
a chaotic region which does not show additional struc- 
tures such as small regular islands or partial barriers. 
Our theory gives an analytical expression for the qual- 
ity factor which is in very good agreement with the full 
numerical simulations of Maxwell's equations. We would 
like to emphasize that our theory can predict quality fac- 
tors also in the regime of large wave numbers where nu- 
merical methods fail due to the exponential increase in 
the quality factor. 
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FIG. 4: (Color online) Inverse quality factor 1/Q for the 
annular microcavity. Shown is the theoretical prediction 
(solid line) which is the sum of the direct tunneling contribu- 
tion (dotted line) and the dynamical tunneling contribution 
(dashed line), and numerical data (dots) for to = 7, ... , 21, 
n = 1 at d = 0.33 and n e g = 2.3. The inset shows a Poincare 
section of the classical phase space, where the critical line 
p = p c is marked (dashed line) . 
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APPENDIX: THE CIRCULAR MICROCAVITY 

For completeness, let us consider a circular microcav- 
ity of radius a in more detail. Inside the cavity the re- 
fractive index is denoted by n\ and outside by 712 < n\ 
(see Fig. [He)). In the classical ray picture trajectories 
stay inside the cavity if their angle of incidence with the 
boundary is larger than the angle of total internal reflec- 
tion arcsin(ri2/?T.i ) • The dynamics is completely regular. 
The circular cavity in TM-polarization is described by 
the Schrodinger equation in polar coordinates (r, if) 

- W 2 ijj{r, ip) = n(r) 2 k 2 i[>(r, tp) (16) 

where n(r) changes from n\ inside to outside the cav- 
ity. The radial and the angular part can be separated, 
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using the ansatz ip(r,ip) — u(r)<f>(ip). We immediately 
obtain <f>(tp) — e lmip , where m e Z denotes the angular 
quantum number. The radial part 



£l ■ J.§;.) ^ VMr)u{r) = h* ll (r) (17, 



describes the motion of a particle in an effective potential 

,2 



V eS {r) = k 2 [l-n(r) 2 } + 



m 



r 2 ■ 



(18) 



where the incoming wave is described by the Hankel func- 

(2,) 

tion of the second kind Hm (n 2 kr) and the scattered 
one is described by the Hankel function of the first kind 
S m Hm\n 2 kr) with a certain scattering amplitude S m . 
A m describes the amount of probability entering the cav- 
ity. Using that the radial solutions and their derivative 
have to be continuous at r — a and that the scattering 
matrix has a pole at a resonance position, this complex 
resonance position k = Rc(fc) + ilm(fc) can be found by 
numerically solving 



Metastable states inside this potential exist for 
m/{n\a) < k < m/(n 2 a) and correspond to states with 
evanescent leakage that, in the ray picture, are fully con- 
fined by total internal reflection. The solutions of the 
radial equation are given as 



u(r) 



A m J m (nikr), r < a 

H$(n 2 kr) + S m H^(n 2 kr), r > a 



(19) 



niJ m+1 (nika)H^ (n 2 ka) = n 2 J m (nika)H { ^ +1 (n 2 ka) 

(20) 

for complex k. From this Qdir = — Re(fc)/[2Im(fc)] is 
obtained. 
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